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RESEARCHES IN THE LUNAR THEORY. 



By Gr. W. Hill, Nyack Turnpike, N. Y. 
Chapter II. 

(Continued from p. 147.) 

The method of employing numerical values, from the outset, in the 
equations of condition, determining the a^, is far less laborious than the 
literal development of these coefficients in powers of a parameter. For com- 
parison with the results just given, we add the calculation of the coefficients 
by this method. The following table gives the numerical values of the sym- 
bols [J, i], [j] and (J), but the division by the quantity 2 (4f — 1) — 4m + m 2 
has been omitted ; it is easier to perform this once for all at the end of the 
series of operations, than to divide each coefficient separately. Hence it must 
be understood that all the numbers in each department of the table are to be 
divided by the divisor which stands at the head of it. 

Coefficients for aj and a_j . 
Divisor = 5.68314 08148 64695. 



[1] = 


0.00861 47842 96261 


[-1] = - 


0.01178 75756 56865 


(1)=- 


0.00623 66553 18347 


(-1) = - 


0.04941 95042 02516 


[1,-2] = 


13.30665 60411 


[-1,-3] = - 


66.98979 68560 


[1, -1] = 


6 32993 22853 


[-1,-2] = - 


28.01307 31002 


[1,2]=- 


10.71949 01593 


[-1, 1] = ~ 


10.96365 06556 


[1,3]=- 


15.10904 80332 


[-1,2]=- 


38.57409 27816 



Coefficients for a 2 and a_ 2 • 
Divisor = 29.68314 08148 64695. 



2 [2] 

2(2) 

[2, -2] 

[2, -1] 

[2,1] 



00205 43632 76229 
0.02909 07097 39048 

14.97672 37558 
9.32666 40103 

13.00326 82750 49 



[- 



2[-2] = 
2 (—2) = ■ 
-2,-4] = 



[-2,-3]=- 
[-2,-1]=- 



0.01834 79966 76898 
0.07227 35586 23216 
108.69586 45706 
63.00980 48251 
8.67987 25398 



G2 
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[2, 3] = — 50.03961 76194 [—2, 1] = — 3 64352 31954 

[2,4]=— 74.07269 86888 [_2,2] = — 19.61044 21261 

Coefficients for a 3 and a_ 3 . 

Divisor = 69.68314 08149. 

[3]=— 0.00113 35729 26473 [— 3]=— 0.00793 43596 

(3)=— 0.01768 33677 (_3) = _ 0.03207 76506 69434 

[3, —1] = 12.99224 12519 [—3, —4] = — 114.67538 20668 

[3, 1] = — 18.10997 74284 [—3, —2] = — 35.57316 33864 

[3,2]=— 41.33769 10334 [—3,— 1]=— 12.34544 97815 

[3, 4] = — 103.14632 67728 " [—3,1]= 1.46318 59580 





Coefficients for a 4 and a_ 4 . 






Divisor = 


125.68314 08 . 




2[4]=- 


- 0.00428 9733 


2[-4] = - 


- 0.01449 0913 


2(4)=- 


- 0.03864 29156 


2(-4)=- 


- 0.06023 435 


[4, -1] = 


16.82502 987 


[-4,-5] = - 


- 182.50817 069 


[4, 1] = - 


- 22.66333 2 


[-4, -3] = - 


- 79 01980 9 


[4,2]=- 


- 51.16496 6 


[-4,-2]=- 


- 42 51817 5 


[4,3]=- 


- 85.50490 2 


[-4, -1] = - 


- 1617823 8 


[4,5] = - 


-171.69968 135 


[-4, 1] = 


6.01654 053 



Coefficients for a g and a_ s . 
Divisor = 197.68314. 

[5]=— 0.00272 9536 (5)=— 0.02896 299 

[5,1]=— 26.99534 4 [—5, -4] = — 138.68780 

[5,2]=— 60.26133 2 [_5 ; _3] = _ 89.421810 

[5,3]=— 99.79795 8 [_5,_2] = — 49.88518 4 

[5, 4] = — 145.60523 2 [_5 > _i] == _ 20.077912 

Coefficients for a 6 and a_ 6 . 
Divisor = 285.68314 . 
2 [6]=— 0.00622 021 2 (—6)=— 0.05640 548 

[6,1]=— 31.21669 [—6,— 5] = — 214.46646 

[6,2]=— 68.99224 [—6, —4] = — 152.69091 

[6, 3] = — 1 13.32666 [—6, —3] = — 100.35648 
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[6,4]= — 164.21995 [—6,-2]=— 57.46319 

[6,5]= — 221.67212 [—6,-1] = — 24.01103. 

These numbers are arranged for carrying the precision to quantities of 
the 13th order inclusive, and to 15 places of decimals. The quantities [j, i\ 
can be tested by differences, if and the divisor with the negative sign are 
inserted in the proper places in the series of numbers ; for it is evident that 
the second differences should be constant. 

The final results are given below, where, in order that the degree of 
convergence of this process may be appreciated, we have given the value 
arising from the first approximation, and then, separately, the corrections 
arising severally from the second and third approximations. It must be borne 
in mind that each of these terms is the numerical value, not of an infinite 
series, but of a rational function of m, and, consequently admits of being 
computed exact to the last decimal place employed, and, in fact, is here so 
computed. Hence any error there may be in these values of the a 4 arises only 
from the neglect of the terms of the following approximations, which, in half 
the number of cases, are of the 14th order, and, in the other half, of the 16th 
order. It is safe to affirm that these cannot, in any case, exceed two units 
in the 15th decimal. 

1st apx , term of 2d order, + 0.00151 58491 71593 — 0.00869 58084 99634 

2d " " 6th " —0.0000001416 98831 +0.0000000615 51932 

3d " " 10th " +0.0000000000 06801 —0.00000 0000013838 

^ = + 0.0015157074 79563, a_ t = — 00869 57469 61540, 

a a 

a 2 . a_2 . 

1st apx.,term of 4th order, + 0.00000 58793 35016 + 0.00000 01636 69405 
2d " " 8th " —0.00000 00006 78490 +0.000000000121088 
3d " " 12th " +0.00000 0000000052 —0.000000000000007 

a, = + 0.00000 58786 56578, a__ 2 = + 0.00000 01637 90486, 

a a 

a 3 . a_ 3 . 

1st apx.,termof 6th order, + 0.00000 00300 35759 + 0.00000 00024 60338 

2d " " 10th " —0.000000000004128 +0.00000 0000000055 

a 3 = +0.0000000300 31632, a_ 3 = +0.0000000024,60393, 
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-4 



&4 • 

1st apx., term of 8th order, + 0.00000 00001 75296 + 0.00000 00000 12284 

2d " " 12th " —0.00000 00000 00028 00000 00000 00000 



a, = + 0.00000 00001 75268, a_ 4 = + 0.00000 00000 12284, 

Of the 10th order, s^= + 0.00000 00000 01107, a^ = + 0.00000 00000 00064, 

a a 

Of the 12th order, a« = + 0.00000 00000 00007, a_ 6 = + 0.00000 00000 00000. 
a a<) 

These give the following numerical expression for the coordinates, 
r cos v = a [1 — 0.00718 00394 81977 cos 2* 
+ 0.00000 60424 47064 cos 4r 
+ 0.00000 00324 92024 cos 6r 
+ 0.00000 00001 87552 cos 8* 
+ 0.00000 00000 01171 cos 10* 
+ 0.00000 00000 00008 cos 12r] , 
r sin v = a [ 0.01021 14544 41102 sin 2r 
+ 00000 57148 66093 sin 4r 
+ 0.00000 00275 71239 sin 6r 
+ 0.00000 00001 62985 sin 8? 
+ 0.00000 00000 01042 sin 10* 
+ 0.00000 00000 00007 sin 12*] . 
On comparison of these values with those obtained from the series in lit, the 
differences are found to be only some units in the 11th decimal. 

The coefficients tend to diminish with some regularity as we advance 
towards higher orders. This is shown by the following scheme of the logar- 
ithms arid their differences : 





A 


A 2 




A 


A 2 


n 97.8561 






98.0091 


- 3.2521 




94.7812 


-2.2691 




94.7570 


2.3165 


-j-9356 


92.5118 


2.2387 


+ 307 


92.4405 


2.2294 


871 


90.2731 


2.2046 


341 


90.2111 


2.1931 


363 


88.0685 


2.1809 


237 


88.0180 


2.1730 


201 


85.8876 






85.8450 
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For verification the following equations were computed, 

2«. [(2» + 1 + m) 2 + 2m 2 ] a,. QS 4 . aj 2 = 1.17141 84591 84518 a* , 
%. (— 1)* (2i + 1) (2i+ 1 + m) a,. [2*. (— 1)* aj 2 = 1.17141 84591 84513 a 2 . . 

The small difference between the numbers is explained by the fact that, in 
these formulae, the quantities a* are, when i is somewhat large, multiplied by 
large numbers ; as, for instance, a 6 by 169. From the average of these two 
results we get 

a = 0.99909 31419 75298 ["£~| * . 

In the investigations of succeeding chapters, the function — plays an 

important part. Hence we will here derive its development as a periodic 
function of r by the method of special values. By dividing the quadrant, 
with reference to t, into 6 equal parts, we obtain the advantage that the sines 
or cosines of the multiples of 2r are either rational or involve V3. -The 

special values of the coordinates and of — , thence deduced, are 



t. — cos V. 

a 

0° 0.99282 60356 45842 

15 0.99378 49245 37167 

30 0.99640 6926450272 

45 0.99999 '39577 40480 

60 1.00358 70309 15127 

75 1.00622 11177 22330 

90 1.00718 60496 23406 



— sin v. 
a 

0.00000 00000 00000 

0.00511 07041 52675 

0.00884 83280 32746 

0.01021 14268 70906 

0.00883 84298 76613 

0.00510 08054 31947 



1.19699 57017 23421 
1.19348 68051 03032 
1.18399 66676 76716 
1.17125 64904 33157 
1.15876 77987 29687 
1.14978 07679 95764 
1.14652 34925 50570 . 



0.00000 00000 00000 
From the numbers of the last column, by the known process, we deduce 

% - - 1.17150 80211 79225 

+ 0.02523 36924 97860 cos 2* 
+ 0.00025 15533 50012 cos 4r 
+ 0.00000 24118 79799 cos 6t 
+ 0.00000 00226 05851 cos 8r 
+ 0.00000 00002 08750 cos 10t 
+ 0.00000 00000 01908 cos 12* 
+ 0.00000 00000 00017 cos 14r . 
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The last coefficient has been added from induction, after which it becomes 
necessary, as is plain, to subtract an equal quantity from the coefficient of 
cos 10r.. Writing the logarithms^ as in the former case, we have, the last 
logarithm being supplied from estimation, 

A A 2 A 3 



98.4020 


-2.0014 






96.4006 




-168 






2.0182 




4-68 


94.3824 




100 






2.0282 




36 


92.3542 




64 






2.0346 




20 


90.3196 




44 






2.0390 




10 


88.2806 


2.0424 


34 





86.2382 
It will be noticed how much slower this series converges than those for the 
coordinates. 

Any information regarding the motion of satellites having long periods 
of revolution about their primaries will doubtless be welcome, as the series 
given by previous investigators are inadequate for showing anything in this 
direction. Hence this chapter will be terminated by a table of the more 
salient properties of the class of satellites having the radius vector at a mini- 
mum in syzygies and at a maximum in quadratures. For this end I have 
selected, besides the earth's moon, taken for the sake of comparison, the 
moons of 10, 9, 8, ... , 3 lunations in the periods of their primaries, and also 
what may be called the moon of maximum lunation, as, of the class of satel- 
lites under discussion, exhibiting the complete round of phases, it has the 
longest lunation. 

In order that the table may be readily applicable to satellites accom- 
panying any planet, the canonical linear and temporal units, that is those 
for which ^ and n' are both unity, will be used. 

From the foregoing methods we obtain : 

For m = — ; 

r cos v ■ = a [1 — 0.011230 cos 2* r sin v = a [ 0.016102 sin 2<r 

+ 0.000015 cos 4t] , + 0.000014 sin 4r] , 

log a = 9.3051648. 
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For m zz 



r cos v — a [1 — 0.014044 cos 2t 

+ 0.0000247 cos 4t] , 



1 

9"' 




sin v - 


- a [ 0.020232 sin 2t 
+ 0.0000230 sin 4r] , 



log a = 9.3326467 



For m zz 



r cos v zz 



a [1 — 0.018061 cos 2t 
+ 0.0000421 cos 4t 
+ 0.00000057 cos 6r] , 



8 ' 

r sin u = a [ 0.026172 sin 2t 
+ 0.0000388 sin 4* 
+ 0.00000048 sin 6r] , 



log a = 9.3630019 



r cos v zz 



For m 

a [1 — 0.02407886 cos 2t 
+ 0.00007760 cos 4? 
+ 0.00000141 cos 6r 
+ 0.000000025 cos 8r] , 



1 

V 



r sin u — 



a [ 0.03516059 sin 2r 
+ 0.00007063 sin 4r 
+ 0.00000118 sin 6r 
+ 0.000000022 sin 8r] 



log a = 9.3969048 



For m zz — 



r cos v zz 



a [1 _ 0.03368245 cos 2t 
+ 0.00015943 cos 4r 
+ 0.000004077 cos 6t 
+ 0.000000097 cos 8r] 



r sin t> zz 



a [ 0.04968194 sin 2t 
+ 0.00014312 sin 4r 
+ 0.000003393 sin 6r 
+ 0.000000084 sin 8r] , 



log a 



9 4352928 . 



For m zz — - ; 
5 



r cos v zz a [1 — 0.05038803 cos 2r 
+ 0.00038127 cos 4t 
+ 0.000014686 cos 6r 
+ 0.000000505 cos 8r] , 



r sm v zz 



a [ 0.07536021 sin 2r 
+ 0.00033582 sin 4t 
+ 0.000012168 sin 6r 
+ 0.000000438 sin 8r] , 



log a zz 9.4795445 . 
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For m 

r cos v = a [1 — 0.08331972 cos 2<r 
4- 0.00114564 cos 4t 
+ 0.00007409 cos 6* 
+ 0.00000404 cos 8*] , 



1 




4 ; 




sin v = 


= a[ 0.12709553 sin 2r 




+ 0.00098090 sin 4r 




+ 0.00006099 sin 6r 




+ 0.00000342 sin 8r] 



log a = 9.5318013. 



For m = 



a [ 0.2542740 sin 2? 
+ 0.0039840 sin 4* 
+ 0.00049306 sin 6r 
+ 0.000070196 sin 8r 
+ 0.000010611 sin 10r 
+ 0.0000016902 sin 12t] , 



r cos u = a [1 — 0.1622330 cos 2t r sin u = 
+ 0.0048920 cos 4r 
+ 0.00059858 cos 6t 
+ 0.000081198 cos 8r 
+ 0.000011873 cos lOr 
+ 0.000001849 cos 12r] , 

log a = 9.5955815. 

For moons of much longer lunations the methods hitherto used are not 
practicable, and, in consequence, we resort to mechanical quadratures. Here 
we shall have two cases. The satellite may be started at right angles to and 
from a point on the line of syzygies, and the motion traced across the first 
quadrant ; or it may be started at right angles to and from a point on the 
line of quadratures, and the motion traced across the second quadrant ; the 
prime object being to discover what value of the initial velocity will make 
the satellite intersect perpendicularly the axis at the farther side of the 
quadrant. 

The differential equations 



£-*i+[?- 8 ]'= ' 



^l 4- 2 — 4- 1. 

dt 2 ^ dt "*" r 3 



= 0, 



give, as expressions of the values of the coordinates, in the first case, 



»=*/!<*—) *-M. s ? 



dt 2 
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and, in the second case, 

*= - 2 X <*-»)* -MX?-*]** • 

y = *, — 2 f'mB— f'f '^dt*. 

J a J oJ o r 

Here the subscript ( ) denotes values which belong to the beginning of 
motion, and ( t ) will hereafter be used to denote those which belong to the end. 

Let v be the velocity, and a the angle, the direction of motion, relative 
to the rotating axes, makes with the moving line of syzygies. In the first 
case then a = 90°, and we wish to ascertain what value of % will make 
a t = 180°. Generally, for small values of % , a x will come out but little less 
than 270° ; but, as v augments, a t will be found to diminish, and, if x does 
not exceed a certain limit, a value of v can be found which will make 
(Tj = 180°. In the second case, in like manner, we seek what value of v 
will make a x = 270°. 

Mechanical quadratures performed with axes of coordinates having no 
rotation possess some advantages, as, in this case, the velocities are not 
present in the expressions of the second differentials of the coordinates. 

Let X and Y denote the coordinates of the moon in this system, and X 
its longitude measured from the line of the last syzygy, from which t is also 
counted. Then the potential function is 

n = — — ^ r 2 -\- -- (X eos t + Ysmt) 2 . 
r 2 2 ' 



F = d dX = - [?" + X ] X+3r cos (;i ~ t] cos *' 



And 

d?X 

dt 

<PY 

it 2 " dY~~ L.r 

Therefore, if we compute p and 6 from 

j> cos d = — — — 2r\ cos (h — t) , 
p sin 6 = — — -f- r sin (a, — t) , 

cPX. 
dt 2 ' 
d 2 Y 



we shall have 



df =iJcos(0 + #), 



df =P - (0 + 0. 



64 
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The needed values of v and <r can be derived from the equations 

v cos (p + t) = — + Y, 
at 

dY 
v sin (a 4- t) — — — — X. 
v y eft 

The developments of the coordinates in ascending powers of t, t being 

counted from any desired epoch, can often be employed with advantage. 

Differentiating the differential equations n times we have 



d n + 2 x „d n + 1 y , d n x d n , 8 N 

— 2 " 4- 3 — — — (r~ 3 x) 

dt H + 2 dt n + 1 ^ dt n dt n v ' 

d n + hi d n + h' d n . 3 N 

^ — — 2 — — (r~ y) . 

dt n + 2 dt n + 1 dt nK J) 



Also 

<!?_ (r -* T \ _ r -3 ^ , _ rf (r~ 3 ) rf'-'a: , n (n — 1) d 2 (r~ s ) d n ~ 2 x 
dt nV '~ dt" ^ dt dt"- 1 ^' 1.2 dfe a ^"- 2 " r " 

with a similar formula for the differential coefficients of r~~' s y . 
The differential coefficients of r -3 , as far as the 4th, are 

dt \ dt^ J dt) ' 

- dF - = ~ 3r (* 5? +y di + W + dt 2 ) + 15r v^ + ^) < 
~7zF~ ~ — \ W + y W + d* 7ft 2 + di dt 2 ) 

, on 7 /" dx . dt/\ ( d 2 x . d 2 y , dx 2 . dy\ 



V dt J dt) ' 



.J4 ( r -:f) 



ft 4 



_ 5 T d 4 # , d 4 ?/ . , <&p <Z 3 # . > dw d s y . fd 2 x\ 2 Q (d 2 y\ 2 ~\ 

- s * l x w + yw+ 4 Ttit> +4 /tw +s w) + s (df)l 

ir 7 / ck , tfy\ / d 3 # , e/ 3 ?/ , Q <&r <Z 2 # , rZv d 2 y\ 
4- 30r- 7 (* (PX 4- « ^ 4- ^ 4- ^Y 

- 5 *-(^+4)*(4f+4?+^+f) 
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By means of these formulae x and y can be expanded in series of ascending 
powers of t, as far as the term involving f\ provided we know the values of 

x, y, -- and -?- corresponding to t ■= . Taking t sufficiently small to make the 
(tv dv 

terms, involving higher powers of t than the sixth, insignificant, as, for instance, 

(1 V (J?/ 

t = 0.05 or t = 0.1 , we can ascertain the values of x, y, — - and — at the end 

' "" dt dt 

of this time. With these values we can again construct new series for x and 
y in powers of t, in which the latter variable is counted from the end of the 
previous time. By repetitions of this process the integration can be carried 
as far as desired. Jacobi's integral, which has not been put to use in the 
preceding formulae, can be employed as a check. 

In case the body starts from, and at right angles to, either axis, the 
coefficients of every other power of t in the series for the coordinates vanish. 
Thus when the axis in question is that of x, the series for the coordinates 
have the forms 

x - x, + A 2 t 2 + A,t A + A 6 t 6 + A 8 t* + . . . , 
y = v t + A s t 3 + A 5 t 5 + A 7 f + A 9 t» + ... . 

By substitution of these values in the differential equations and the equating 
of each resulting coefficient to zero we arrive at the following equations ; 

1 . 2 A 2 = 2% + 3x — x^ 2 , 

2.3A, = — 4A 2 — aff\ , 

3 . 4 A, = 6J 3 + 3A 2 + \ x^ (3vl + 4x A 2 ) , 

4 . 5 A 5 = — 8A 4 + ~ a%\ {v\ + 2x A 2 ) — x^ s A 3 , 
5.6A, = 10A 5 + 3A 4 + I a*T* (6»„^ 8 + 4x A 4 + 3A\) 



2 

v ^o" 6 («o + 2x A 2 ) (5v 2 + Qx A 2 ) , 



Q IK 

6 . 7 A 7 = — 12A 6 + | xv\ (2v A, + 2x A, + A$) — { x^X {< + 2x Q A 2 f 

+ - 2 - x^ (vl + 2x A 2 ) A 3 — xi 3 A 5 , 
7.8A 8 = UA 7 + 3A % + | xr 4 (6^ 5 + 4x A 6 + 6^ 2 ^ 4 + ^§ 

— 4-arjT 8 K + 2.r .4 2 ) (10iv4, + 8x ^ 4 + ^|) 
+ i 5 6 aS" 8 K + 2 *oA 2 y (7vl + 8^ 2 ) 
+ | x - 5 (2« J 3 + 2avi 4 + At) A t , 
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8 . 9 A, = - 16A S + 4 a^\ (2v A s + 2x A 6 + 2A,A, + Al) 

15 



j AV r $ (% + 2;r J 2 ) (2v A s + 2x A, + ^) 



85 



+ ii rt K + 2ff ^) 8 + 2 av 6 (2v A 3 + 2^ 4 + Jf) A 3 
— j W 1 (< + 2x A 2 ) 2 A 3 + ]- xv 5 (vl + 2d-„^ 2 ) J 5 — x^A 7 . 
By means of these relations each A can be derived from all the A which 
precede it. 

When the axis is that of y, the series have the forms 

x = v t + A,f + A 5 t 5 + A 7 t 7 + A t° + ..., 
y = y,+ A 2 f- + A^ + A 6 f + A s t s + . . . . 
And the equations, determining the coefficients A, are 
1 . 2 J 2 = — 2v — y^ 2 , 
2.3A 3 = 4A 2 + 3v — y^v , 

3 . 4 ^ 4 = — QA, + j #r* (3» 2 + % J 2 ) , 

4 . 5 J 5 = 8 A, + 3A 3 + ~ yr % (vl + 2y A 2 ) — y^A s . 

The equations are not written as far as in the former case, as it is evident 
they may be derived from the preceding group by putting y in the place of 
x , reversing the signs of the first terms, and removing the term 3A n _ 2 from 
the equations, which give the values of the A of even subscripts, into those 
which give the values of the A of odd subscripts, after having augmented the 
subscript by unity. 

The velocity of the moon of maximum lunation vanishes in quadratures, 
and when v = the preceding seines become, putting y^ 3 = a , 

' = * [- t at> + (as • " m "*) e + (~ 2626 ■ + as «' + m «') ' 

( 1 1 2 1 „ 47 A s 
+ Vi81440 a ~" 12096 a * + 45360 a ° + 9072 a V *" 

/ 1 317 2 _ 13 3 _ 1.0403 4,_947_ A 1 

+ Vl9958400 a+ 9979200 a 120960 a 4989600 a + 237600 a / ' J' 

/ 1 „ 1 2 , 13 3 73 



1 V10080 1008 + 1120 a 5040 /* 
1 _17_ 2 _ _1_ 3 4603 
907200 a + 90720 a 756 a + 453600 ~ 400 



+ ( 1 a + ^L a 2 - 1 a 3 + 46 ° 3 a 4 - 3 a 5 ") H 

+ V 907200 + 90720 756 a+ 453600 400 / J " 
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These series suffice for computing the values of x and y with the desired 
exactitude when t is less than 0.3. 

This special case of the moon of maximum lunation will now be treated. 
As there seems to be no ready method of getting even a roughly approximate 
value of y , we are reduced to making a series of guesses. I first took y = 0.82 ; 
tracing the path to its intersection with the axis of x , a t , which ought to be 
270°, came out 261° 29' 47".9. A second trial was made with y = 0.7937 ; 
the result was a, = 267° 37' 8".3. Again a third trial with y = 0.7835 gave 
(tj = 269° 41' 13".3. The principal data acquired in the three trials are given 
in the following lines : 

nl T n. l "'l/l - Maximum 

9 dt fit X Variation. 

0.8200 0.972430 —0.339523 —0.288149 —1.927275 261° 29' 47".9 44° 57' 4" 
0.7937 0.908207 —0.290945 —0.089184 —2.144832 267 37 8 .3 46 39 36 
0.7835 0.884782 —0.274324 —0.012170 —2.227928 269 41 13 .3 47 17 21. 

1 denotes the time employed in crossing the quadrant, and the last column 
contains the maximum value of the angular deviation of the body from its 
mean direction as seen from the origin, that is, the direction it would have 
had, had it moved across the quadrant with a uniform angular velocity about 
the origin. 

A check may be had on the accuracy of the computations by mechanical 
quadratures. We determine the value of the constant 2C which completes 
Jacobi's integral from the coordinates and velocities, both at the beginning 
and at the end of the motion, for each of the three trials. The result is 

Second value. 

2.43901 

2.51987 
2.55261 . 

We can now apply Lagrange's general interpolation formula to these 
data, and, regarding a x as the independent variable, inquire what are the 
values which correspond to a t = 270°. The numbers of the first trial must 
be multiplied by + 0.014861 ; those of the second by — 0.210190 ; those of 
the third by + 1.195329 , and the sums taken. The results are 



y . 


First value. 


0.8200 


2.34902 


0.7937 


2.51985 


0.7835 


2.55265 



dx t dy x 

dt ' dt 

0.781898 0.881160 0.271798 —0.000083 —2.24093 2.55788 47° 23' 12". 



a, T «, ®X 1 Cty i cy n Maximum 

dt dt Variation. 
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fit* 

That — — does not rigorously vanish is due to the employment of only three 

terms in the interpolation ; for the same reason the value of 2C does not quite 

dy ± 
dt 



agree with that obtained from the values of x t and — . To make all these 



elements accordant we add 0.00009 to the value of -{— . 

dt 

A table of approximate values of x and y, derived roughly from the data 
afforded by the process of mechanical quadratures is appended: they will 
serve for plotting the orbit. 



t. 


X . 


y- 


t. 


X . 


y> 


t. 


x . 


y ■ 


0.00 


— .0000 


+ .7819 


0.30 


— .0148 


+ .7080 


0.60 


— .1177 


+ .4748 


0.02 


.0000 


.7816 


0.32 


.0180 


.6978 


0.62 


.1294 


.4519 


0.04 


.0000 


.7806 


0.34 


.0215 


.6869 


0.64 


.1418 


.4277 


0.06 


.0001 


.7790 


0.36 


.0256 


.6752 


0.66 


.1547 


.4022 


0.08 


.0003 


.7767 


0.38 


.0301 


.6629 


0.68 


.1680 


.3752 


0.10 


.0005 


.7737 


0.40 


.0351 


.6499 


0.70 


.1818 


.3466 


0.12 


.0009 


.7701 


0.42 


.0407 


.6361 


0.72 


.1956 


.3162 


0.14 


.0015 


.7659 


0.44 


.0468 


.6216 


0.74 


.2095 


.2839 


0.16 


.0022 


.7610 


0.46 


.0534 


.6063 


0.76 


.2230 


.2496 


0.18 


.0032 


.7554 


0.48 


.0607 


.5902 


0.78 


.2359 


.2131 


0.20 


.0044 


.7492 


0.50 


.0686 


.5733 


0.80 


.2475 


.1745 


0.22 


.0058 


.7432 


0.52 


.0771 


.5555 


0.82 


.2575 


.1339 


0.24 


.0076 


.7347 


0.54 


.0863 


.5369 


0.84 


.2653 


.0913 


0.26 


,0096 


.7265 


0.56 


.0961 


.5172 


0.86 


.2704 


.0474 


0.28 


.0120 


.7176 


0.58 


.1066 


.4965 


0.88 


.2718 


.0027 



The following is the table of the numerical values of the quantities of 
principal interest belonging to the moons mentioned at the beginning of 
this paragraph. In the first line stands the earth's moon, having very 
approximately 12-rA- lunations in the period of its primary. In the last line 
is the moon of maximum lunation. The quantities belonging to the moon of 
two lunations have been somewhat rudely inferred from the numbers in the 
adjacent lines. 
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N umber of 

Lunations 

in period 

of Primary. 


Radius 
Vector 

in 
Syz3 T gies. 


Radius 
Vector 
in Quad- 
ratures. 


Ratio. 


Velocity 

in 
Syzygies. 


Velocity 
in Quad- 
ratures. 


Ratio. 


20. 


1 

m 


r - 


r x . 


7\_ 


V . 


«!. 






tO 59 


0.17610 


0.17864 


1.01446 


2.22295 


2.16484 


0.97386 


6.50888 


10 


0.19965 


0.20418 


1.02271 


2.06163 


1.97693 


0.95892 


5.88686 


9 


0.21209 


0.21813 


1.02849 


1.98730 


1.88501 


0.94853 


5.61562 


8 


0.22652 


0.23485 


1.03678 


1.90904 


1.78250 


0.93372 


5.33873 


7 


0.24342 


0.25543 


1.04934 


1.82721 


1.66572 


0.91162 


5.05535 


6 


0.26332 


0.28167 


1.06969 


1.74333 


1.52851 


0.87677 


4.76409 


5 


0.28660 


0.31699 


1.10605 


1.66247 


1.35953 


0.81777 


4.46103 


4 


0.31232 


0.36897 


1.18138 


1.60111 


1.13480 


0.70876 


4.13277 


3 


0.33235 


0.45973 


1.38329 


1.62141 


0.79387 


0.48962 


3.72018 


2 


0.302 


0.684 


2.26 


2.00 


0.18 


0.09 


2.89 


1.78265 


0.27180 


0.78190 


2.87676 


2.24102 


0.00000 


0.00000 


2.55788 



Maximum 
Variation. 



0°35' 6" 

55 21 

1 9 33 

1 29 58 

2 53 
2 50 49 
4 18 37 
7 17 

14 34 14 
37 21 

47 23 12 



In regard to this table we may notice the following points. The moon of 
the last line is the most remarkable : it is, of the class of satellites considered 
in this chapter, (viz., those which have the radius vector at a minimum in 
syzygies, and at a maximum in quadratures,) that which, having the longest 
lunation, is still able to appear at all angles with the sun, and thus undergo 
all possible phases. Whether this class of satellites is properly to be pro- 
longed beyond this moon, can only be decided by further employment of 
mechanical quadratures. But it is at least certain that the orbits, if they do 
exist, do not intersect the line of quadratures, and that the moons describing 
them would make oscillations to and fro, never departing as much as 90° 
from the point of conjunction or of opposition. 

This moon is also remarkable for becoming stationary with respect to the 
sun when in quadrature ; and its angular motion near this point is so nearly 
equal to that of the sun that, for about one-third of its lunation, it is within 1° 
of quadrature. From the data of the table we learn that such a moon, circu- 
lating about the earth, would make a lunation in 204.896 days. 

We notice that the radius vector in syzygies of this class of satellites 
arrives at a maximum before we reach the moon of maximum lunation. This 
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maximum value is very nearly, if not exactly, i, when measured in terms of 
our linear unit, and thus is a little less than double the radius vector of the 
earth's moon. It occurs in the case of the moon which has about 2.8 luna- 
tions in the period of its primary. 

The radius vector in quadratures augments continuously as the length of 
the lunation increases, as also does the ratio of these radii, until, in the moon 
of maximum lunation, the radius in quadratures is but little less than three 
times that in syzygies. 

The velocity in syzygies does not continuously diminish, but attains a 
minimum somewhere about the moon of four lunations, and afterwards aug- 
ments so that, for the moon of maximum lunation, it does not differ greatly 
from the velocity of the earth's moon in syzygies. On the other hand the 
velocity in quadratures constantly diminishes. 

The maximum value of the variation augments rapidly with increase in 
the length of lunation, so that, in the moon of maximum lunation, it exceeds 
an octant, or is more than 80 times the value which belongs to the earth's 
moon. 

In the adjoining figure are constructed graphically the paths of the 
earth's moon, of the moons of four and three lunations, and of the moon of 
maximum lunation. The moons in the first lines of the table have paths 
which approach the ellipse quite closely, but the paths of the moons of the 
last lines exhibit considerable deviation from this curve, while the orbit of 
the moon of maximum lunation has sharp cusps at the points of quadrature. 

( To be continued. ) 



Plate 111. 
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